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We study the giant magnetoresistance of nanogranular magnets in the presence of an external 
magnetic field and finite temperature. We show that the magnetization of arrays of nanogranular 
magnets has hysteretic behaviour at low temperatures leading to a double peak in the magne- 
toresistance which coalesces at high temperatures into a single peak. We numerically calculate the 
magnetization of magnetic domains and the motion of domain walls in this system using a combined 
mean-field approach and a model for an elastic membrane moving in a random medium, respectively. 
From the obtained results, we calculate the electric resistivity as a function of magnetic field and 
temperature. Our findings show excellent agreement with various experimental data. 
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INTRODUCTION 

Recent experimental studies of magnetic nanoarrays [l], 
i, S i, S, S 0, i S, 113, EH, [II El, El revealed the giant 
magnetoresistance (MR) effect. At low temperatures the 
hysteresis of the magnetization was found to lead to a 
splitting of the resistance peak 0, 0, H, H, E^ whereas 
at temperatures above ^ lOOK the magnetization hys- 
teresis disappears and the MR shows a single peak at 
zero external magnetic field. The interest in these sys- 
tems is motivated both, by the important technological 
promise [ll| and by the opportunity of applying the ideas 
developed for granular materials for the description of 
doped manganite systems [l6. 17|. 

In this paper we investigate the giant MR of magnetic 
nanoarrays and develop a model for the field and temper- 
ature dependent magnetization hysteresis and the sam- 
ple resistivity. We consider temperatures below the Curie 
temperature of the individual grains, T-!. For T <T^ the 
nanoarray can be either in a superferromagnetic (SFM) 
state with the magnetic moments of the grains aligned by 
ferromagnetic interactions and the systems has a mul- 
tidomain structure [11], or, at higher temperatures, in 
a superparamagnetic (SPM) state, where magnetic mo- 
ments of individual grains (superspins) are mutually dis- 
oriented by thermal fluctuations. In real systems the 
ferromagnetic coupling can appear due to a competi- 
tion of exchange and dipole-dipole interactions. The for- 
mer is always present which favors ferromagnetic order, 
whereas the latter alone would lead to an antiferromag- 
netic order or a striped phase depending on the orienta- 
tion of the (super-) spins. Recently a detailed analysis of 
nano-granular arrays of magnetic particles revealed that 
short-ranged exchange interaction stabilize the superfer- 
romagnet and the ferromagentic interaction is provided 
by "glue particles", which are basically remains of the 
sample preparation process of nearly atomic size [3]. 

The electronic transport within a single domain in the 
SFM state or in the SPM state is controlled by the aver- 



age mutual alignment of the superspins [l^ [see Eq. ([T]) 
below] . In addition, in the SFM state the contribution of 
the domain wall dynamics in an applied magnetic field 
and the associated change of the size of domains with 
different bulk magnetizations on the MR has to be taken 
into account. 

Magnetoresistance in granular materials was stud- 
ied theoretically in Refs. [H, 21, 
ever, these considerations were classical, 



before. How- 
using a spin- 
dependent scattering approach with a phenomenological 
spin-diffusion length, which is defined as the distance a 
spin-polarized conduction electron travels before it un- 
dergoes a spin-flip collision, for regular arrays of monodis- 
perse spherical grains. In contrast to these previous con- 
siderations, we take Coulomb blockade effects, multiple 
electron co-tunneling, and domain wall dynamics into ac- 
count in the present paper . Moreover, our approach is 
applicable for arrays with position disorder, size distribu- 
tion of the grains, quenched mesoscopic disorder due to 
impurities within the grains as well as in the embedding 
matrix. On top of the imperfections of the system, also 
annealed disorder in the form of thermal fluctuations of 
the magnetic moments is included. 



ELECTRON TRANSPORT 

We consider two- and three-dimensional arrays of con- 
ducting magnetic nano-size granules, weakly coupled via 
tunnel junctions. Granularity gives rise to Coulomb 
blockade effects, controlling the low temperature conduc- 
tivity at T <C i^c {Ec = e? /{na) is the charging energy 
of a single granule, with e being the electron charge, n 
is the sample dielectric constant, and a is the granule 
size. It can be as large as several hundred Kelvins which 
is exactly the experimentally relevant temperature inter- 
val) [sl] . Then the conductivity of a nanogranular array 
composed of superspins in a single super-domain in the 



2 



hT 



domain t 



domain i 




(2) 



FIG. 1: Illustration of an idealistic array of ferromagnetic 
nano-particles with two magnetic super-domains in the SFM 
state. Our model also includes effects of quenched disorder 
induced by e.g. size and position fluctuations of the grains. 
The domain wall is considered to be thin compared to the 
domain size; the bulk-magnetizations m-f ,| of the two domains 
have opposite signs. The bold wavy line shall only illustrate 
the model of an elastic string which represents the real domain 
wall. In general domain walls are several grains thick where 
the average orientation of the superspins changes gradually 
from T to J,. The dynamics of the domain wall is described by 
an elastic object driven by the external magnetic field /ioH 
moving in the z-direction. 



SFM (or in the SPM) state is given by 

a{T, rr?) ~ g°(l + H^m^) exp(- V^o(m2)/T), (1) 

where is the tunneling conductance of the SPM state 
with no magnetic field, 11 the polarization factor of the 
grains, m? — {cos 6) the normaUzcd bulk magnetization 
{9 is the angle between the magnetic moments of adja- 
cent superspins and (. . .) stand for averaging over a single 
domain), and To{m'^) = To[l - (Co/a) ln(l -hH^m^)] with 
To = e^/(K^o) and ^0 being the characteristic energy scale 
and the inelastic localization length, respectively. Equa- 
tion ([T|) describes the contribution from the co-tunneling 
process es g overning hopping conductivity in granular ma- 
terials [23[ which holds as long as the electron hopping 
distance exceeds the size of a single grain, a, i.e. as long 
as T < Tq. 

The MR in nanogranular magnets (see Fig. [T|) stems 
from (i) the field dependent bulk magnetization and (ii) 
the domain wall motion in the SFM state. Normally the 
thickness of the domain walls is negligible as compared to 
the domain size itself, allowing to model their dynamics 
under the external magnetic field by an elastic membrane 
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driven through a random medium \24 
We also exclude domain nucleation processes within a 
domain involved and assume that the domain walls are 
sufficiently spatially separated or the external magnetic 
field changes fast enough in order to avoid effects of do- 
main wall interaction (for more details, see discussion 
at the end of the paper ). Thus, only the change of 



domain sizes associated with the domain walls dynam- 
ics contribute to the electron transport but not to the 
electron scattering at the walls [s^]. Incorporating the 
field dependent bulk magnetization and domain dynam- 
ics into the description of electron transport of Ref . [31 , 
we construct a theory of the giant magnetoresistance in 
magnetic nanoarrays. 

MAGNETODYNAMICS 

First, we discuss the bulk magnetization in the SPM 
state and/or in a single domain in the SFM state, using 
the mean-field (MF) equation [i^ 

m{h, T) = L{a) , with a = [rjh + 3L(a)]/T , (2) 

where L(a) = coth(a) — l/a is the Langevin function, 
and we introduce the following dimcnsionless parameters: 
h = H/Ho, T = T/T^, and rj = ^iHo/ikBT^). Here 
Hq is a characteristic magnetic field strength and /i is 
a material dependent density of magnetic moments. In 
experiment [6| the typical values are iiqHq « 0.2T, « 
lOOK, and r/ « 1. Equation ^ is vahd in the SPM state, 
and in the SFM state for h > hc{T) with hc{T) < being 
the critical field where the MF approach breaks down. At 
this field a{h) in Eq. ^ becomes discontinuous. 

Next, we turn to the SFM state. The equation of mo- 
tion for an over-damped D-dimensional interface profile 
z(x, t), describing the domain- wall, embedded in a d = 
D + 1 dimensional system has the following form [3^ [sOl : 

~ = TV^z + fit) + .g(x, z) + 7?(x, t) . (3) 
7 ot 

Here 7 and P denote the mobility and the stiffness con- 
stant of the interface, respectively, and f(t) is the time 
dependant driving force (/(i) = 2fiofJ,BH(t) with /is 
and ^oH(t) being the Bohr magneton and magnetic 
field, respectively). Random forces ^(x, z) are defined by 
(g(x, z)) = and (5(x, z)g(x', z')) = 5^(x - x')Ao(|z - 
z'l). We further assume that Ao(|z|) is a monotonically 
decreasing function of z for z > which decays to zero 
over a finite distance £, [s^ . Thermal fluctuations are de- 
scribed by the random noise term 77(x, t) with {t]{:x., t)) = 
and (r?(x,i)77(x',i')> = 2^ 5^(x - x') (5(t - i')- For the 
details of model ^ we refer the reader to Ref. [2^ . 

If the system is driven adiabatically, i.e. if it is always 
in a steady state even if the driving force is changing, 
the mean velocity v(t) = {dz{:x.,t)/dt)^ [i^ shows a de- 
pinning transition at zero temperature or creep at finite 
temperatures 25[ . As we are interested in the finite tem- 
perature regime and aim to describe a finite system, we 
restrict the domain wall by |z(x, f)| < Lz, and define the 
magnetization parallel to the magnetic field as 



MJT,h) 



[(1 + z)m-^ + (1 — z)mi] 



(4a) 



3 



T=0.01 




T<1 



T>1 



FIG. 2: Magnetization hysteresis, Eq. (|4b|l . for different tem- 
peratures below T = T /T" < 1 and above the Curie temper- 
ature of the sample according to Eq. ([2]). The arrows show 
the direction of hysteresis with changing magnetic field. The 
dimensionless magnetic field h is given in units of hoHq. 



where Mg is the saturation magnetization, z = ^(T, h) = 
{z{x,t))^ /Lz the dimensionless mean displacement of 
the domain wall, and m|,x = m|,x(T, /i) the bulk mag- 
netization for the two domains with opposite mean di- 
rection with mx(T, ft,) — — m|(T, — ft,). The normalized 
sample magnetization is therefore 



m,(T,ft)EEA/||(T,ft)/A/,, 
i.e. ms(T, ft) e [-1, 1]. 



(4b) 



SIMULATIONS 

We solve Eq. ^ numerically for various temperatures 
and for f{t) slowly changing between -/max and /max, 
such that non-adiabatic effects like the hysteresis of the 
domain- wall velocity v{t) can be neglected [13, S]- We 
introduce dimensionless coordinates and time in the fol- 
lowing way: The units of (simulation) time tq and space 
Ao are chosen such that 7rTo = Aq and that the dimen- 
sionless random pinning forces g ~ Toj/X^g are uniformly 
distributed in [—1/2,1/2]. The dimensionless driving 
force is ft = /Ao/(to7) and the dimensionless temper- 
ature t = To7/A^+^r. 

We now discretize Eq. ^ in x-directions in lat- 
tice sites (for D = I the lattice Laplacian is given by 



Zi+i + Zi-i — 2zi, i = 1, . . . , N.). The stochas- 



tic forces g{x., z) are constructed by choosing a random 
number in [—1/2, 1/2] at z-positions apart by i. In be- 
tween the forces are linearly interpolated which results 
in a Gaussian g-g correlator with variance £, [3l[. 

The dimensionless sample magnetization ms(T, ft) for 
different temperatures T and ftmax = 3 is shown in Fig.[2l 
By comparison with experimental data of Ref. 0] the 
dimensionless simulation temperature T corresponds to 
~ T/10 and ft = 1 to a magnetic field of 0.2T. Notice, 
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FIG. 3: Magnetic field dependent magnetoresistance for vari- 
ous temperatures. Temperatures are given in units of T". At 
low temperatures T < 1 the magnetization hysteresis leads to 
the double peak structure of the MR (left y-axis). In the SPM 
regime (T > 1) the double peak is coalesced into a single more 
pronounced peak (right y-axis). The resistivity is normalized 
to the maximal resistivity. 



that for \h{f)\ > ftc(T) « 2.5 at T = 0.01 the contri- 
bution of the opposite sign domain with respect to ft 
can be neglected for the sample magnetization. Since in 
real experiments the external magnetic field ft cannot be 
changed adiabatically/infinitly slow, the driving force in 
the simulation is also changed very slowly but still oscil- 
lating, resulting in hysteretic behaviour of ms(T, ft) even 
at finite temperatures. 

Using the results for the bulk magnetization and the 
domain wall dynamics below we can now calculate 
the magnetoresistance for our model. In the SFM regime 
we have to distinguish the two cases of parallel and per- 
pendicular resistivity with respect to the domain wall: 
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(1 + z)a{J, m2) + (1 - z)a(T, mj) '^^^^ 



We note that p^(T,G) = p||(T,0) since m|(T, 0) = 
TO^(T, 0). In the SPM state the domain structure does 
not exist. Thus, the electron transport is governed by the 
behaviour of bulk magnetization described by Eq. ^ . 

Using typical values for the parameters 11^ = 0.3, 
^o/a = 1, Tq/T^ = 10, we can calculate the magnetic 
field dependent magnetoresistance, p{T, ft) = [p_l(T, ft) -f- 
P||(T,ft)]/2 in the SFM state and p(T, ft) = l/cr(r,m2) 
in the SPM state where m(T, ft) is derived from Eq. 
The resulting curves for various temperatures are shown 
in Fig. [3l At low temperatures, T < 1, a clear double 
peak in the resistivity of the system can be seen which 
is a direct result of the magnetization hysteresis. The 
observed peak depends on the sign of the external mag- 
netic field change, i.e. the right (left) peak for increas- 
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ing (decreasing) field. In the SFM state the resistivity 
can only be calculated for \h\ < hc{T), where the bulk 
magnetization varies only slightly such that there is no 
visible difi^erence between p± and pn . For T > 1 the MF 
approach holds for all h. 

The maximum change of the resistivity with magnetic 
field can be quantified using the expression Ap(T) = 
(/5max(T) - /3oo(T))/pinax(T), where Pmax(T) is the maxi- 
mum (peak) resistivity and poo (T) the " saturation" resis- 
tivity for the fully magnetized state [m^(T, /i — > ±00) = 
1]. The factor by which p changes is given by 7 — 
Ap(T)/[l-Ap(T)]. 

In the SFM state we get Ap{T = 0.01) w 0.032 or 7 = 
3% and Ap{T = 0.1) w 0.107 or 7 = 12%. The decrease 
of the effect towards lower temperatures is due to the 
fact that the bulk magnetization is less field dependent. 
In the SFM state we can write an explicit expression 
for Ap{T), using Eq. (P), given by Ap(T) = 1 - (1 + 
n2)-ig-ViVT+V'7'o(i)/T_ Pqj, t = 2 we get 7 = 78%. 
Note, that within the field range of Fig. [3] the curve for 
T = 2 does not saturate. 

We now compare our results with reported experimen- 
tal data. In Refs. [l|, 0, H, 0, H, @, the magnetization 
curves for magnetic nanoarrays were measured mostly at 
high temperatures - showing a weak hysteresis effect. In 
Ref. the magnetization data for weakly coupled 10- 
nm Co-nanocrystals has been obtained at 5K where the 
magnetization hysteresis has a width of approximately 
O.IT leading to the double peak structure of the resis- 
tivity. The resulting maximal resistivity change 7oxp is 
about 6%. The double peak structure in the resistiv- 
ity was observed in several experiments @, 0, 0, 0| , with 
Ap in the range 2% to 10%, which is in good agreement 
with our results. However, mostly the single peak in the 
giant magnetorestance in the SPM state at high temper- 
atures (from lOOK to room temperature) has been stud- 
ied 0,li,i,S,H[iil,[i3l. It would be interesting to com- 
pare our predictions for the magnetoresistance p{T , h) to 
more precise measurements at low temperatures near the 
sample Curie temperature. 



DISCUSSION 

There are several characteristic length scales in the 
present problem. Some of them are relevant for the 
description of the magnetization m^(T,/i) while others 
for the electron conductivity aij'^m?). To describe the 
magnetization behaviour in the SFM state we used the 
elastic model for domain walls driven in a random en- 
vironment, Eq. ([3]): comparing the elastic (first) and 
random force (third) terms in the r.h.s. of Eq. ([3]) one 
can show that the Imry-Ma or Larkin length scale is 
given by ip « [(ff )VAo(0)]i/('*--°). On this scale the 
weak random forces accumulate to a value comparable to 
the elastic force. On smaller length scales the domain- 



wall is essentially flat and above the interface becomes 
rough. Therefore the domain size Ldom in nanogranu- 
lar magnets should be larger than Lp in order to use 
the model for domain wall dynamics. Thus, our con- 
sideration is valid for the following sequence of length 
scales, Ldom > Lp > a. On the other hand the rele- 
vant length scale for electron conductivity is the hop- 
ping length rhop(r) ~ a^/To/T. Since the energy scale 
Tq ^ Ec > T the hopping length cannot be smaller than 
the size of a single grain, rhop(T') > a- At the same time 
the relation between rhop (T) and idom is arbitrary. 

The mapping of simulation temperature and magnetic 
field slightly depends on the rate of change of the exter- 
nal driving force h{t) in Eq. ([3]) - if it is changed more 
slowly, hysteresis effects will become smaller. If the ex- 
ternal driving force h{t) is changed faster and the do- 
main wall motion becomes non-adiabatic - in the sense 
that e.g. the velocity hysteresis becomes pronounced pTj 
" several effects should be mentioned [28|, |32| : First the 
magnetization hysteresis gets wider since the domain wall 
cannot follow the changes of the magnetic field anymore, 
leading to a greater separation of the double peaks of 
p{T,h). Furthermore the change in the magnetization 
from m = ±ltom = =pl will become smeared out lead- 
ing to a broadening of the resistivity peaks. In addition 
to these visual effects a phase shift between the magnetic 
field and magnetization appears. The shape of the mag- 
netization curves near the saturation points depends also 
on boundary effects - in the simulation z was simply re- 
stricted by the system border, but no additional surface 
pinning effects were taken into account. The qualitative 
picture remains, however, the same. 

In our model we consider only one domain wall. There- 
fore, we comment under which conditions this model is 
justified, in particular why dipolar domain wall interac- 
tion can be neglected. Since the dipolar interaction be- 
tween walls is of Coulomb type it is proportional to l/£^ 
(in 2D), where £ is the typical distance between domain 
walls [33[ . From this we can introduce a typical time scale 
for the domain wall interaction as r = M/q'^^ where q is 
the charge of the walls. If the frequency of the external 
magnetic field is larger than 1/r the contribution of do- 
main wall interaction can be neglected. In other words 
if domain wall fluctuations due to interaction are slow 
(e.g. for large spatial distances between them) compared 
to the motion induced by the driving force, the former 
motion has no influence on the sample resistivity. We 
assume this in our model and therefore consider only one 
domain wall. The spiked domain structure observed in 
LMOKE micrographs in Ref. [3] cannot be described in 
detail by our model. However, the effect of these struc- 
tures on the magnetoresistance can be assumed to be 
small such that a coarse grained picture with only one 
wall will describe the transport correctly, at least on a 
qualitative level. We remark, that samples with a single 
domain wall can be created, see e.g. Ref. [35| where an 
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ultrathin Pt/Co/Pt film was prepared with one moving 
wall, and model Q was successfully applied to describe 
the magnetic susceptibility of nanogranular magnets [28j . 

In general, also dipolar interactions within one domain 
wall have to be considered in the Hamiltonian for the 
wall, leading to an additional term in Eq. ^ . This term 
describes the interaction between dipoles via Coulomb 
forces. Since the dipolar axis is usually parallel to the 
wall this additional term represents only a renormaliza- 
tion of the elastic term which can be seen in Eq. (25) of 
Ref. [s^]. Therefore its influence is hidden in our phe- 
nomenological constant F in Eq. 

In conclusion, we have studied the giant magnetore- 
sistance of nanogranular magnets in the presence of an 
external magnetic field and finite temperature. We have 
shown that the magnetization of arrays of nanogranular 
magnets has hysteretic behaviour at low temperatures 
leading to a double peak in the magnetoresistance which 
coalesces in the SPM state to a single one. We have 
numerically calculated the magnetization of this system 
based on the model for domain walls moving in a random 
medium combined with a MF approach for the bulk mag- 
netization. Fig. [21 Using the simulation results for the 
magnetization we have calculate the electric resistivity 
p{T, h) as a function of magnetic field and temperature 
in the Ohmic regime, Fig. [3l Our findings are in a good 
agreement with published experimental data 0, S, H 0| • 
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